Second-order effects on the hyperfine structure of P states of alkali-metal atoms 
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We analyze second-order Ml-Ml and M1-E2 effects to the hyperfine structure (HFS) of the 
lowest energy P states of alkali-metal atoms arising from the coupling of the two (J = 1/2, 3/2) 
fine-structure levels through the hyperfine interaction. We find these effects to be especially sizable 
in Li, leading to a 9a correction to the most accurate reported experimental value of the A(P 1 / 2 ) 
HFS constant of 7 Li [D. Das and V. Natarajan, J. Phys. B 41, 035001 (2008)]. For the remaining 
alkali-metal systems, the results tabulated within may be referenced as higher precision is sought in 
experimental determination of the HFS constants. 

PACS numbers: 32.10.Fn,31.15.A- 



(N 



I , 



cd : 

o ■ 

r/3 ■ 



> 

O 

O 
oo 
O 



% 



I. INTRODUCTION 

As experimental accuracy improves, the interpretation 
of measurements may require refined theoretical analy- 
sis. Measurements of the hyperfine structure (HFS) is 
one of such examples. The HFS of atomic systems arises 
from the coupling between the atomic electrons and the 
nuclear spin. Generally speaking, in the approximation 
that J, representing the electronic angular momentum, 
remains a "good" quantum number, the HFS can be ac- 
curately described by the conventional (first-order) HFS 
constants A,B,C, . . . [1;]. Respectively, these constants 
describe the electronic interaction with the nuclear mag- 
netic dipole (Ml), electric quadrupole (E2), magnetic 
octupole (M3),. . . ; these interactions are collectively re- 
ferred to as the hyperfine interaction (HFI). Likewise, 
with J being a "good" quantum number, experimental 
measurements of the HFS intervals, along with known 
values of the nuclear moments, can be used to definitively 
determine these HFS constants. However, for atomic 
states that are part of a fine-structure manifold, the 
nearby fine-structure levels may provide sizable contam- 
ination to the J-purity [2,y, 0, 0- In these situations, to 
make the connection between HFS interval measurements 
and HFS constants, it becomes necessary to consider the 
HFI coupling between fine-structure states [lj. In terms 
of perturbation theory, this arises in the second-order in 
the HFI and is dominated by the Ml-Ml and M1-E2 
effects. 

In this paper, we provide results for calculations of the 
second-order corrections to the HFS for the lowest en- 
ergy P states of the naturally occurring alkali-metal iso- 
topes. These isotopes, along with their respective nuclear 
properties, are displayed in Table [IJ The HFS constants 
A(Pi/ 2 ), A(P 3 / 2 ), and B(P 3 / 2 ) have been experimentally 
measured for all of these isotopes (see Refs. @, la 0| and 
references within). In a recent paper, Das and Natara- 
jan [a claim a significant improvement in experimental 
uncertainty over earlier works for most of these constants. 

Due to its relatively small fine-structure splitting, the 
second-order effects of the HFI are most noticeable in 
Li, the lightest of the alkali-metal atoms. Over thirty 
years ago, in experimental determination of the HFS con- 



TABLE I: Nuclear properties of stable alkali-metal isotopes. 
I v represents the nuclear spin and parity. Nuclear dipole 
and quadrupole values are given in terms of nuclear magne- 
tons (un) and barns (b), respectively. All data is taken from 
Ref. @1. 



Isotope 


r 


M (jJ>N) 


Q(b) 


°Li 


i+ 


0.8220473(6) 


-0.00082(2) 


7 Li 


3/2" 


3.2564625(4) 


-0.0406(8) 


23 Na 


3/2+ 


2.2176556(6) 


0.1045(10) 


39 K 


3/2+ 


0.39150731(12) 


0.0585(6) 


41 K 


3/2+ 


0.21489274(12) 


0.0711(7) 


85 Rb 


5/2~ 


1.3533515(8) 


0.277(1) 


87 Rb 


3/2- 


2.751818(2) 


0.134(1) 


133 Cs 


7/2+ 


2.5829128(15) 


-0.00355(4) 



stants, Orth et al. [2, 3] recognized the need to consider 
the mixing of the fine-structure P states of Li. How- 
ever, more recent measurements of the 6,7 Li A(P 1 / 2 ) 
constants [a Q have overlooked this important effect. 
Namely, we find that the second-order correction to the 
7 Li A(P 1 / 2 ) constant causes a shift which is a full order 
of magnitude more than claimed uncertainty of Das and 
Natarajan J6J. For the remaining alkali- metal isotopes, 
we find that the second-order effects are not significant 
at the current levels of experimental precision for the 
HFS constants. However as higher precision is sought, 
it may become necessary to include these effects. The 
compilation of second-order corrections provided herein 
may provide useful complementary data for future high- 
precision experiments on these isotopes. 

This paper is organized as follows. In Section [Til we 
present pertinent equations for the HFI and correspond- 
ing first- and second-order energies obtained from a stan- 
dard perturbation theory analysis. In Section ITTT1 we pro- 
vide a reformulated tensorial analysis of the perturbation 
theory; this formalism provides more insight into the ro- 
tational symmetries of the various second-order contri- 
butions. We compile HFS expressions for the P states 
of the alkali-metal atoms in Section HVl In Section [V] 
we tabulate our results, followed by concluding remarks 
in Section IVIl Furthermore, we include an Appendix 



containing notations and expressions for spherical tensor 
operators appearing in Sections [II] and IIIII 



II. THE HYPERFINE INTERACTION 



over scalar products of spherical tensors: 



TT \ * rpn rpe 

#HFI - 2^ 1 k ■ J-ki 



(1) 



In this Section, we recapitulate the basic properties 
of the hyperfine interaction and the application of per- 
turbation theory to analyze its effects on atomic struc- 
ture. This Section is an abbreviated review of the pre- 
sentation provided in Ref. [4j. We employ the notation 
(Tk)i to represent the expectation value of the zero- 
component operator of spherical tensor (of rank k) Tk 
in the "stretched" state \I,Mi = I). This is related to 
the reduced matrix element through the expression 



(T k )j = 



I k I 
-I I 



(i\\n\\i). 



The hyperfine interaction can be expressed as a sum 



where T^ and T| are spherical tensors of rank k (k > 0) 
acting in the nuclear and electronic spaces, respectively. 
As i^HFi is a scalar operator in the combined space, it is 
convenient to work in the conventional basis formed by 
coupling nuclear, \IMj), and electronic, \jJMj), states, 

\ 1 IJFM F )= Y, C^ JMj \IM I )\^JM J ), (2) 

with 7 encapsulating remaining electronic quantum num- 
bers and the coupling coefficients being the Clebsch- 
Gordon coefficients. A matrix element of the HFI be- 
tween these basis states is given by 



(<y'IJ'F'M' F \H H Fi\7lJFM F ) = SftS. 



M'Mf 



(-1) 



I+J' + F 



E 



j'il\ {wmunnhJ)- 



First-order and second-order energy corrections to the state described by electronic quantum numbers 7 and J are 



W; 



(1) 



(ryIJFM F \H HFI \~fIJFM F ) = (-1 



J+ 



I. ^ 



I J F 

k 



(7||T,"||/)(7J||T||| 7 J), 



(3) 



t'j' 

7 ' j' 



{'jIJFM F \H H Fi\iIJ'FM F ){iIJ'FM F \H HFI \- ( IJFM F ) 



(-1) 



j-j' 



E, 



iJ 



E. 



itE 



yj 



k±,k 2 



I J F 
J' I fci 



Ey'J' 



I J F 

J' I k 2 



ai|r fc " 1 ||/)(/||T fc " 2 ||/)< 7 J||r fc e i ||7'J')<7'J'IKIl7J). (4) 



r 



Here the summations exclude the case (7V) = (7J); 
this will be implicit in similar summations to follow. 
We note here that, due to large energies associated 
with nuclear excitations, we treat the nuclear states as 
"good" quantum states. The reduced matrix elements 
(IllT^HJ) appearing here are associated with the nuclear 
moments. Specifically, the magnetic dipole (k = 1), elec- 
tric quadrupole (fc = 2), and magnetic octupole (fc = 3) 
are given by 



cally, through the C constant, these are given by 



A = -h(T?)i(TZ)j = ±ni{T?,j, 
B = MT^tiT^j = 2Q(Ti)j, 

c = {t^^tDj = -n{T§)j, 



(5) 



Q = 2(T 2 ") 7 , 

n = -(r 3 ")/. 



Furthermore the product of reduced matrix elements 
(I\\Tg\\I){-yJ\\T%\\~fJ) appearing in Eq. © correspond to 
the conventional first-order hyperfine constants. Specifi- 



For electronic states within a fine-structure manifold, the 
leading second-order effects are due to mixing with the 
nearby states within the manifold (i.e., 7' = 7). Further- 
more, these effects are dominated by the dipole-dipole 
(Ml-M 1) and dipole-quadrupole (M1-E2) terms. The 
constants r) and £ have been used in Refs. [J, [f| to pa- 



rameterize these effects and are given by 



=F 



T 



(/ + i)(2/ + i) Ai 2 |( 7 j||r 1 e || 7 J±i)| 2 



^7 J 



(J+1)(2J+1) 21 + 3 



^ 7 j±i 



J V 21 - 1 

M Q( T ^||Tf||' T J±l>(' T J||T 2 e || 7 J±l) 



F. 



7 J 



F. 



'7J±1 



> (6) 



(7) 



(If both J ± 1 levels exist, it would be necessary to dis- 
tinguish two independent r/'s and £'s.) Explicit formulas 
for matrix elements of the electronic tensors T£ are given 
in Ref. d. 



III. TENSORIAL ANALYSIS OF 
SECOND-ORDER EFFECTS 

In this section we reformulate the second-order con- 
tributions to the energy. We classify the second-order 
terms according to the underlying rotational symmetry 
in the nuclear and electronic spaces. As a result, spe- 
cific second-order terms are connected to the first-order 
constants A, B, C, . . . It becomes clear from a physical 
standpoint why, for instance, the second-order constant 
•q cannot affect experimental determination of the C con- 
stant, as proven by brute manner in Ref. [4(. However, 
the formalism here is also farther reaching, as it can eas- 
ily be extended to other second-order terms past r\ and 
to higher orders. 

We begin with an effective Hamiltonian, H*j, which 
produces the exact hyperfine energies when acting 
on the (unperturbed) coupled basis states, Eq. @, 
i.e. Hfj\jIJFM F ) = W F \-iUFM F ). We decompose 
this effective Hamiltonian into contributions of increas- 



H. 



(2) 



~lJ 



ing orders of the HFI, H$ = H^j 

with associated energy contributions given by W F ' 

(-/IJFM F \H {n j ) \-/IJFM F ). From Eq. ©, we may infer 



(m) 



(1) 



that the first-order Hamiltonian takes the form H^j = 
i^HFi- Furthermore, from Eq. (J4j) , we infer the second- 
order Hamiltonian to be 

H 1 J = -ffHFI-RyJ-ffHFI = ^2 ( T ki ' T ki) R lJ ( T k 2 ' T k 2 ) ' 

fci,fc 2 

(8) 
where the latter expression is obtained by using Eq. JTJ) to 
represent Hhfi- The operator R 7 j here is the resolvent 
operator; it acts in the electronic space and is given by 
the expressions 



R 






E 

yj'F'Mt, 



\iLJ'F'M' F )(iIJ'F'M' F \ 



E. 



~iJ 



E. 



yj> 



It is important to realize that the resolvent operator be- 
haves as a scalar operator under rotations in the elec- 
tronic space. 

We may recouple the tensor operators of the second- 
order Hamiltonian, Eq. ([8]), to isolate parts acting in the 
nuclear and electronic spaces. The resulting expression 
is 

\k\+k 2 +k 



H. 



(2) 



~lJ 



E 

k\ ,k2,k 



(-1)* 



x{T£®T fe " 2 } fc .{T^i? 7 /r fe e 2 } A 



I k L-fci "-* "T J ^2 Jfc • 

(See the Appendix for notational descriptions; see 
Ref. [lCW ] for general relations involving spherical tensors.) 
The second-order correction to the energy is given by the 

diagonal matrix elements of H^j in the coupled basis: 



IF, 



(2) 



= (7lJFM F \H%}\'yIJFM P 



(2) 



= (-1 



J+J+F 



\fcl+fc 2 +fc 



(9) 



E (-D 

x<J||{T£®T£}J|J) 

^{lJ\\{T c kl ®R^T c k2 } k \\lJ)- 

Comparing this to Eq. ([3]), we can see clearly that all 
/c-dependent terms enter the first- and second-order en- 
ergy expressions in an identical manner. However, be- 
cause -S/hfi does not contain a monopole contribution, 
the index k in Eq. ([3]) is limited by k > 0, whereas the 
summation here is inclusive of the case k — 0. Analyzing 
the F-dependent factors entering this expression, namely 
the phase factor and six-j symbol, under the condition 
of k = gives 



•1) 



I+J+F 



I J F 
J I 



v /(27+l)(2J + l)' 



and we see that F-dependence is lost. Thus, the second- 
order k = terms here only cause an overall shift to the 
HFS energy levels, and do not affect the interval spacing 
between levels. 

If all second- and higher-order effects are negligible, 
then the HFS intervals are determined completely by 
the terms (/||T fe ™||/)(7J||T fc e ||7J) appearing in Eq. Q; 
we mention again that these terms are related directly 
(see Eq. {[5])) to the conventional hyperfine constants A 



(k = 1), S (jfc = 2), C (k = 3), 



As TT 



and T^ are 



spherical tensors, the specific value of k represents the 
underlying rotational symmetry in the nuclear and elec- 
tronic spaces, respectively. Experimental measurements 
of the HFS interval spacings may be used to determine 
these constants. However, if the second-order effects can- 
not be neglected, then the terms (7||T^||I)(7J||r|||7J) 
can no longer be determined by the intervals. The "con- 
stants" which could be determined would also depend on 
second-order effects: 



(7||T fc "||/)(7J||T fc e ||7J)+^(-l) 



ki-\-k 2 -\-k 



fcl,/L"2 



x(/||{r &1 ®T£} fc ||/)<7J||{T fc e i ® RyjT^lfrJ). 



We note that second-order dipole-dipole (fci = k 2 = 1) 
terms only affect determination of the A and B constants, 
as two rank-one tensors can only be coupled to form a 
tensor of rank k < 2. This agrees with the proof given 
in Ref . ,4] , in which the C constant was shown not to be 
affected by the second-order constant r\. However, here 
we extend this conclusion to dipole-dipole terms which 
mix states outside of the fine-structure manifold as well. 
Similar conclusions can also be drawn here; for exam- 
ple, second-order dipole-octupole terms do not affect the 
determination of the A constant, as a rank-one and a 
rank-three tensor cannot be coupled to form a tensor of 
rank k = 1. 

The reduced matrix elements appearing here can be 
simplified, giving 



in terms of the HFS intervals 5W ( F J) = W { p ] - W { pl t and 
the constants r\ and £. 

For 6 Li (/ = 1) we have 



A(Pi/ 2 ) 
MP3/2) 
B(P 3/2 ) 



-SW {Pl > 



36^ 



12V3 



= Vifr } + lw^ ] + 1 



6^3/2 



3"" 3/2 



10 5 /2 
5*5/2 



1 



72 



36 ' 



120^ 
1 



c. 



20\/3 



c 



{I\\{Tl®Tl} k \\I) = {-lY i+k V2kTl 



k\ k 2 k 
III 



{I\\Tt\\I){I\\Tl\\I), 



frJII^ ® RyjT&^J) = (-l) 2J+k V2kTT 



E 

7'./' 



k\ k 2 k 
J J J 1 



( 7 J||T fc e i || 7 'J')<7'J , ||T fc e 2 ||7J) 



E, 



iJ 



E, 



~f>j> 



These expressions may be used in Eq. ([9]) , at which point 
the summation over k may be carried out analytically by 
using a well-known six-j sum rule; the resulting expres- 
sion is identical to Eq. ([J}. 

The tensorial analysis of second-order HFI effects pre- 
sented in this Section could be applied to higher-order 
effects as well, with similar insights following. For exam- 
ple, we could find that experimental determination of the 
D constant is not affected by third-order dipole-dipole- 
dipole terms, due to the fact that three rank-one tensors 
cannot be coupled to form a composite tensor with the 
rotational symmetry of a rank-four tensor. 



For the isotopes 7 Li, 23 Na, 39 K, 41 K, and 87 Rb (I 
3/2) we have 



MP1/2) 
A(P V2 ) 

B(P 3/2 ) 



C(P 



3/2J 



Wi*" 5 - - 



2 
] 
20 



c, 



90 15V5 



25 



100 



C, 



180 lSOv^ 

4 1 5 2 20 3 

1 1 „ 

30 ' 20V5 
±8W[ P ^ ] -LsW^ + i^' 1 



80 
400 V5 



100 



400 



IV. EXPRESSIONS FOR P-STATES OF ALKALI 
METAL ATOMS 

In this Section we concern ourselves with the HFS 
equations of the lowest energy P-states of alkali-metal 
atoms, namely the isotopes of Table |U Specifically, 
we are considering the nP\i 2 and nP 3 i 2 states, with 
n = 2,3,4,5,6 for isotopes of Li, Na, K, Rb, and Cs, 
respectively; the specific n will be implicit hereafter. 
As inferred from the arguments of the previous Section, 
to determine the (first-order) constants A,B,C,... by 
measurement of HFS intervals, it is necessary to have 
a knowledge of the higher-order effects. Here we as- 
sume that these higher-order effects are sufficiently de- 
scribed by the second-order dipole-dipole and dipole- 
quadrupole constants r\ and £. These are defined in 
Eqs. (|6I7P ; one convenience of these definitions is that 
77 and C are identical for both the P\/ 2 and P 3 / 2 states 
(i.e., r](Pi/ 2 ) = 77(^3/2) = V an d similar for £). Below we 
compile specific expressions for the first-order constants 



For b5 Rb (/ = 5/2) we have 



MP1/2) 
MPs/2) 

S(P 3/2 ) 



C(P 



3/2; 



1 



sw. 



(P1/2) 



3 
3 

50 2 

1 



315^ 105v/30 



c, 



swy 3/ 



(^3/2) 



64 
525' 



JWo 



70 



-8W t 



(P3/2) 



630 



V 



525V30 



c 



1 m,CV.) 8 sTirlfya) , 15 



-SW. 



1 



().--; 



'/- 



14V30 



21 

c, 



SWi 



-SW, 



(P3/2) 






h s <"" 



28 



+ 1 >-f"> 



280V30 



Finally, for 133 Cs (I = 7/2) we have 



A(P 3/2 ) = -5W 3 

1 



(^3/2) 



Iswpri + ^Vf 3/2) 



21 



120 



1 



VI 



tC, 



1512 ' 1260 

b(p 3/2 ) = -\sw:,~ : --<^ 



8 



120^ 



1_ 1 

"36 ??+ 120 C ' 



C(P 3/2 ) = ±8W< P >'' ) -±Swi P *» ) + -L6Wf''> ) 



32 

+ 48CT 



24 



480 



We note that in all of the above cases, the r\ (£) con- 
tribution to A(P 3 / 2 ) is suppressed by a factor of 1/2 
(—1/10) compared to its contribution to A(P l / 2 ). 



V. RESULTS 

In order to calculate the second-order constants r\ and 
£ from Eqs. (|6I7[) . we must first generate the P\/% and 
P 3 / 2 electronic states. This involves solving the electron 
correlation problem. We employ an ab initio relativis- 
tic coupled-cluster method which includes single, dou- 
ble, and triple excitations from the lowest-order Dirac- 
Hartree-Fock state. While both core and valence sin- 
gle and double excitations are included, the triple ex- 
citations involve simultaneous excitation of the valence 
electron with two core electrons. In other words, the 
triple core excitations are not incorporated in the many- 
body wavefunction. We refer to this method as the 
CCSDvT method. The most sophisticated approxima- 
tion within the CCSDvT method is described in our 
previous work [ll| for three-electron Li. For heavier 
alkalis (beyond Li) the approximation of Ref. [Il[ be- 
comes computationally expensive and we only keep the 
lowest-order terms on the right-hand-side of the valence 
triples equation (T V [D V ] and T V [D C ] terms in notation of 
Ref. [HI]). For the electronic reduced matrix elements 
entering Eqs. (|6l7p . this method gives results with accu- 
racies better than 0.1% for the lightest (Li) and better 
than a few percent for the heaviest (Cs) alkali-metal sys- 
tem. 

Table [IT] displays our computed values of the reduced 
matrix elements of the electronic tensors Tf and T 2 be- 
tween P states, along with experimental values for the 
fine structure energy splitting. The variation in these 
values for different isotopes of the same atomic system 
is below the level of precision shown. Combining these 
values with the nuclear data for each isotope, Table [TJ 
we obtain the second-order constants rj and £. These 
are displayed in Table IIII1 Also displayed in this Ta- 
ble are the "corrections" to the HFS constants ^(Pi/2), 



A(P 3 / 2 ), B(P 3/2 ), and C(P 3 / 2 ) deduced from 77 and C- 
Here each correction is regarded as the difference be- 
tween the actual constant and the measured constant 
based only on first-order equations. These can be in- 
ferred from the equations in Section ITVl e.g., for 133 Cs: 
AA(P 1/2 ) = (1/75% + (1/126)C. 



TABLE II: Computed reduced matrix elements of the elec- 
tronic tensors Tf and T£ between P states and experimental 
fine structure intervals AEf s =Ep^ /0 — Ep yl „\ fine structure 
intervals are taken from Ref. 112 



^3/2 



"1/21 



Atom 



(P3/2m\\p 1/2 ) 



{P3/2\\Ti\\P 1/2 ) 

(MHz/b) 



AE (S 
(cm- 1 ) 



Li 


19.0 


11.9 


0.3366 


Na 


12.1 


59.4 


17.1963 


K 


-12.6 


-103 


57.600 


Rb 


20.1 


213 


237.598 


Cs 


29.2 


313 


554.11 



The most accurate reported value of the 7 Li A(P 1 / 2 ) 
constant based on experimental measurement of HFS in- 
tervals is that of Das and Natarajan, Ref. Q, wherein 
they give a value A(P 1/2 ) = 46.024(3) MHz. However, 
this value is based on first-order HFS equations and is 
thus completely negligent of higher-order effects. From 
Table IIII1 we see that the second-order effects would 
cause a sizable shift of AA(P 1 / 2 ) = 27.0 kHz to this 
value. This shift is a full order of magnitude larger than 
the claimed accuracy (3 kHz) of the constant itself. Das 
and Natarajan also claim accurate results for A and B 
constants for the P states of all isotopes in Table UJ with 
the exception of the 6 ' 7 Li P 3 / 2 and 41 K Pi/2,3/2 states. 

Figure [TJ displays various reported values of the 7 Li 
A(P 1 / 2 ) constant and their associated error bars. Along 
with the value from Das and Natarajan, the Figure in- 
cludes experimental values from Orth et al. [3j and Walls 
et al. @ ■ The earlier value of Orth et al. accounted for 
second-order effects, whereas that of Walls et al. did not; 
both claim an accuracy of the same order as our com- 
puted A^4(P 1 / 2 ). Also included in this Figure are the- 
oretical values recently reported by our group [ll| and 
Yerokhin [Tj|]. This Figure also displays "corrected" val- 
ues using our computed AA(P 1 / 2 ) for the applicable cases 
of Das and Natarajan 0] and Walls et al. 0. We note 
that this correction causes these two values to shift far- 
ther away from the values of Orth et al. |3j , Derevianko et 
al. [11| , and Yerokhin [13J ; we do not attempt to explain 
the resulting large discrepancy. 

For the remaining HFS constants of Li, the second- 
order correction again proves to be sizable, though not 
as pronounced as in the 7 Li A{P 1 / 2 ) case. For the 6 Li 
A{P\/2) constant, Das and Natarajan claim an accuracy 
of the same order of magnitude as our predicted second- 
order correction AA(P 1 / 2 ) = 4.01 kHz. 

Orth et al. 0,3 performed HFS measurements of 6:7 Li 
in the presence of strong fields and used a complex fit- 
ting scheme to obtain the HFS constants. In addition 



TABLE III: Second-order constants r\ and £ and corresponding second-order corrections to the HFS constants. All values are 
in kHz. x[y] denotes x x lO^. 



Isotope 



'/ 



c 



AA(P 1/2 



4.01 

2.70[+l] 

2.14[-1] 

6.94[-3] 

4.25[-3] 

4.01 [-2] 

1.08[-1] 

4.11[-3] 



AA(P 3/2 ) 
2.01 

1.41[+1] 

4.03[-2] 

-7.20[-5] 

-2.38[-4] 

-2.34[-3] 

8.32[-3] 

2.38[-3] 



AJ?(P 3/2 ) 
4.02 

8.34[+l] 
3.92[-l] 
7.54[-3] 
3.89[-3] 
6.28[-2] 
1.53[-1] 
9.70[-2] 



Ag(g 



3/2 J 



"Li 

7 Li 

23 Na 

39 K 
41 R 

85 Rb 
87 Rb 
133 Cs 



1.45[+2] 
2. 53 [+3] 
9.25 

9.33[-2] 
2.81[-2] 
8.76 [-1] 
2.87 
3.51 



-2.02[-l] 
-3.41 [+1] 

3.72 

1.98[-1] 

1.32[-1] 

2.68 

2.57 
-6.69[-2] 



-3.82[-2] 
4.16[-3] 
2.21 [-4] 
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FIG. 1: (Color online) Values of the HFS constant A(P 1/2 ) for 
7 Li. The heavy (hollow) squares denote experimental (theo- 
retical) values reported in the respective reference, and the 
heavy diamonds denote the associated corrected values (of 
the immediately preceding value) using second-order correc- 
tion AA(P 1/2 ) = 27.0 kHz from Table Hill when applicable. 
The uncertainty of A J 4(P 1/ / 2 ) is negligible compared to the 
experimental uncertainties in these cases. No uncertainty is 
given for d. a Ref. i; b Ref. [|; c Ref. [l]; d Ref. [H; e Ref. [l|. 



TABLE IV: Computed off-diagonal constant A(P 3/2 ,P 1 / 2 ) 
for 7 Li compared to non-relativistic and experimental values. 
Units are MHz. 

A(P 3 /2,Pi~ 



CCSDvT ab initio; this work 
non-relativ. ab initio; Ref. [14] 
experiment; Ref. [3J 



11.9 

11.87853 

11.823(81) 



r\ = 22.46 kHz was used |17|. In a following paper, 
Ref. [lq . this effect was reevaluated using a more so- 
phisticated third-order approach, resulting in a value of 
i] = 2.09 kHz. Our present value (77 = 3.51 kHz), based 
on the more complete CCSDvT method, is 68% larger 
than the latter value. 

For the heavier alkali-metal atoms past Li, all of the 
predicted second-order corrections (Table Hn|) are one to 
six orders of magnitude smaller than the uncertainty in 
the most precise reported experimental value of the re- 
spective HFS constant. Thus, at the current levels of pre- 
cision, the neglect of second-order effects in these heavier 
systems is acceptable. However, as experimental tech- 
niques improve and higher precision is sought, it may 
become necessary to consider these higher-order effects. 



to the traditional constants, this scheme also yielded an 
off-diagonal constant A(P 3 / 2 , P\/2)- This constant can 
be related to the off-diagonal matrix element of the Tf 
tensor by 



1 



A(P 3/2) P 1/2 ) = - 7 = J(P 3/2 ||7T||P 1/2 



M, 



12 

This constant, along with the fine structure interval, pa- 
rameterizes the mixing between the two fine structure 
levels. In Table IIVI we compare our value of this off- 
diagonal constant for 7 Li (deduced from the matrix el- 
ement of Table |nj with the non-relativistic value from 
Ref. [1J| and experimental value reported by Orth et 
al. [3j. Our (relativistic) value is in close agreement with 
both values. 

The second-order dipolc-dipole correction was consid- 
ered for 133 Cs in Ref. [15J, where a theoretical value of 



VI. DISCUSSION 

In this paper, we have demonstrated for the P states 
of the alkali-metal atoms that the second-order effects 
cannot simply be neglected in deducing HFS constants 
from high-precision measurements of the HFS intervals. 
This is exemplified by the 7 Li A(Pi/ 2 ) constant, where 
the second-order effects are shown to cause a 9cr shift to a 
recently reported value. Furthermore, we have compiled 
values for the leading second-order HFS effects caused by 
the dipole-dipole constants r\ and the dipole-quadrupole 
constants £ for all naturally occurring alkali-metal iso- 
topes. These values can be used in conjunction with 
high-precision measurements of the HFS intervals of the 
alkali- metal atoms to determine the HFS constants. 

This work was supported in part by the National Sci- 
ence Foundation. 



APPENDIX A: COUPLING OF SPHERICAL 
TENSORS: NOTATIONS AND DEFINITIONS 

Here we provide a description for the notations used in 
Sections [TT] and IIIII along with basic definitions involving 
coupling of spherical tensors. For a more complete de- 
scription of spherical tensors, including specific formulae 
for recoupling of spherical tensors and matrix elements 
of coupled tensor operators, see Ref. [10J | . 

A spherical tensor of rank k, P&, is a set of 2k + 1 oper- 
ators. We denote the individual components (operators) 
as Pkq, with the index q taking all integer values from 
— k to k. Two spherical tensors, Pk t and Qk 2 , may be 
coupled together to form a composite spherical tensor of 
rank k. The rank of the composite tensor is limited by 
I hi — &2 1 < k < k\ + &2 • We denote the coupled tensor 
by {Pkt ® Qk 2 }k'i ^s components are given by 



where C 



{Pk^QkAk. 



kq 



/ j ^k 1 q 1 -k 2 q 2 ^k 1 q 1 Qk 2 q 2 



(Al) 



kq 
k 1 q 1 ;k 2 



represents Clebsch- Gordon coefficients. 



The scalar product of two tensors of equal rank is defined 
as 



= {-lfV2k + l{Pk®Qk} 
= £(-l) 9 P feg Q fe - g . 



Pk-Qk 



If one of the tensors in Eq. (JATJ) is a scalar, the compo- 
nents of the coupled operator are simply 

{Po ® Qk} kq = PoaQkq ; 

{Pk <g> Qo} k q = PkqQ0Q ■ 

When P (Q) here is understood to be a scalar, the 
composite tensor is then sufficiently represented by PQk 

(P k Q). 
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